Localized Bosonic Modes in Superconductors 
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We show that a locahzed bosonic mode acts as a new type of "defect" in s- and d^2_y2-wa,ve 
superconductors. The mode induces bound or resonance states, whose spectral signature are peaks 
in the superconductor's density of states (DOS). We study the peaks' shape and energy as a function 
of temperature and the mode's frequency and lifetime. We identify several characteristic signatures 
of the localized mode that qualitatively distinguishes its effects from those of magnetic or non- 
magnetic impurities. 
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The study of local defects or impurities in supercon- 
ductors has attracted significant experimental 0,01 and 
theoretical Q |^ interest over the last few years. 
These studies have proven particularly important for 
further elucidating the nature of the superconducting 
(SC) pairing mechanism in unconventional supercon- 
ductors. In particular, recent scanning tunnelling mi- 
croscopy (STM) experiments provided a detailed picture 
of the frequency and spatial dependence of defect induced 
resonance states in the high-temperature superconductor 
(HTSC) Bi2Sr2CaCu208+5 0, and the triplet supercon- 
ductor Sr2Ru04 0- Several theoretical scenarios for the 
physical origin of these impurity states have been pro- 
posed, ranging from electronic scattering off classical im- 
purities 3, 5] to the onset of Kondo-screening 4,]. 

In this Letter we investigate a new type of "defect" , a 
localized bosonic mode, and study its effect on the local 
electronic structure in s- and dj.2_y2-wave superconduc- 
tors. The study of localized modes, and in particular 
phonon modes in normal metals, has been of great in- 
terest in the context of molecular spectroscopy 0j S 
(or inelastic electron tunnelling). It was recently argued 
that a localized magnetic mode could also arise from un- 
screened magnetic moments in a dj.2_y2-wave supercon- 
ductor Eo]. The nature of the mode is only relevant for 
the present study to the extent that it determines the 
form of its coupling to the electronic degrees of freedom. 
While theoretical studies 0| have so far only focused 
on the Born-limit of weak electronic scattering off the lo- 
calized mode, we show that the strong-scattering, unitary 
limit gives rise to an abundance of novel phenomena. In 
particular, we show that the presence of a localized mode 
leads to the appearance of fermionic bound or resonance 
states whose spectral signature are peaks in the SC DOS. 
We study the peaks' energy and shape as a function of 
the mode's characteristic frequency ujq, lifetime F"^, and 
temperature. It is of particular interest that the peaks 
move to higher energies with increasing or F, but are 
shifted to lower energies with increasing temperature, T. 
In addition, the temperature evolution of the peaks dif- 
fers qualitatively if its frequency is smaller or larger than 
Ao — loq, where Aq is the maximum SC gap. Moreover, 



we show the mode also induces a "dip" in the DOS at 
frequencies ±(Ao-l-a'o)- This dip, together with the tem- 
perature evolution of the DOS, qualitatively distinguishes 
the effects of a localized mode from those of non-magnetic 
or static magnetic impurities These characteristic 

differences therefore provide an important tool for future 
experiments to further study and clarify the nature of 
"defects" . 
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FIG. 1: (a) The full electronic Greens function, (b) Lowest 
order vertex corrections. 

Starting point for our calculations is the Hamiltonian 
H = Xl*^k4^^Ck,^ + ^Ak4^4^-f /i.e. 



k k 

+uJoblbQ + ^ ga (bl + 5o) 

k,q,cj 
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Here, cj., 5q are the fermionic and bosonic creation oper- 
ators, respectively. We consider a two-dimensional elec- 
tronic system with normal state dispersion ek and SC 
gap Ak. The bosonic mode with frequency lvq is local- 
ized at r = 0. We allow for a spin dependence of the 
scattering vertex, g^- , and for simplicity consider only an 
on-site interaction. In the following, we study a single 
bosonic mode with no internal degrees of freedom, i.e., 
spin S — 0. The extension of our results to S > 1, 
in which case electronic spin- flip scattering is allowed, is 
strai ghtf orward and will be discussed in a future publica- 
tion The retarded bosonic Greens function Dr{uj) 
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in frequency space is given by 



where a depends on T and tjQ, and is determined via 

duj 



(2) 



1 + 2nf,(wo) 



■ nb(cjo)Im£'i?(c^) . (3) 



We assume that and the mode's lifetime, F"^, are 
determined intrinsicaUy such as in the case of phonon 
modes 0- Possible feedback effects on cjq and V due 
to the fermion-boson coupling are accounted for (to the 
extent that they do not induce a frequency dependence 
of r) by considering loq and T to be effective tempera- 
ture independent input parameters of our theory. The 
unperturbed (clean) fermionic Greens function is given 
by G^^(k, iujn) = [iw„ro - CkTs] (Tq + Akr2cr2 , where 
and Ti are the Pauli-matriccs in spin and Nambu-space, 
respectively. Changes in the fermionic Greens function 
due to the scattering of the localized mode are accounted 
for by summing the infinite series of diagrams shown in 
Fig. which yields 

+Go{r, 0, LUn)S{LOn)Go{Q, r', w„) . (4) 



Here 



^(wn) = E(w„) i - Go(0, 0,w„)i;(u;„) 



(5) 



and E is the local fermionic self-energy given by 



T.{uJn) = 60(0, 0,W„ - U,n)D{Vrn) Tag . (6) 



with g being determined by the spin dependence of ga- 
We show later that vertex corrections (see Fig. ^) are 
negligible for the range of parameters considered here. 
Based on earlier studies 0, llfl , we expect that the inclu- 
sion of non-crossing diagrams in the calculation of G does 
not qualitatively change our results and therefore reserve 
its discussion for a future publication The DOS, 

N — All + A22 with Aii{r,uj) = — 2ImG,;i(r, + iS) is 
obtained numerically from Eqs.Q-ljni) with 6 — 0.2 meV. 

As a frame of reference, we first consider a localized 
mode in an s-wave superconductor with normal state dis- 
persion ek = fc^/2m — /z, m~^/Ao = 15, Fermi momen- 
tum kp — 7r/2, and a momentum independent SC gap 
Aq. For a spin-independent coupling, g = gl, no bound 
states are induced, similar to the case of a non-magnetic 
impurity in an s-wave superconductor In contrast, if 
the coupling is spin-dependent the mode induces a bound 
state, with poles in the S'-matrix leading to two peaks in 
the DOS. For simplicity we take g-^ = g, and gi — ~g, 
i.e., g = ga^, and in Fig. [3 present the DOS at r = 
for r = 0^ as a function of loq and temperature. For 
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FIG. 2: DOS at the site of the localized mode (r = 0) with 
r = 0^. (a) T = and several values of uo- The curves for 
ujo ~ and ljq = 2Ao are multiplied by 0.1. Inset: arrows 
indicate the "dip" in the DOS. (b) ljo = 0.5Ao and several T. 



T = 0, cjQ = (see Fig. the frequencies ztfli, of the 
bound state peaks are given by rif,/Ao = D-/D+ where 
D± = 1 ± (gm/A)'^. This result is similar to the ener gy o f 
a bound state induced by a static magnetic impurity |l2l| . 
Note that increasing g, i.e., stronger scattering, leads to 
a decrease of ilb- For ujo,T ^ the bound state frequen- 
cies are shifted from fif,, and in the limit T <C Ao are 
determined by a transcendent equation. It is too cum- 
bersome to be presented here, but can be simplified in 
the limit ri{,, <JC Aq, where we find to leading order in 
luq/Aq and nsiuJo) 



r2fc(T,cjo) = fib + (^o/t^ - n_B(a;o)Ao 



(7) 



Note, that the localized mode effectively behaves as an 
"oscillating vector" with Ising symmetry (for a single 
bosonic degree of freedom) and frequency wq- When loq 
increases (and in particular, when it becomes compara- 
ble to the rate of electronic scattering, ^i/h) subsequent 
scattering events begin to cancel each other; the effective 
scattering becomes weaker, and the peaks are shifted to 
higher energies. In particular, for luq —^ 00 the bound 
state vanishes and the DOS becomes that of the unper- 
turbed systems |ll|. This picture provides an explana- 
tion for our analytical results in Eq.Q and fully agrees 
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with the numerically obtained DOS shown in Fig.[2K for 

T = and several values of wq. A non-zero wq also leads 
to a "dip" in the DOS at — Aq +u!q, as shown in the 
inset of Fig. Oi. This dip is a characteristic signature of 
the localized mode, which distinguishes it from a static 
magnetic impurity. The origin of the dip lies in the poles 
of the retarded self-energy, T,ft.{k,fl) at |r2| = Ek + ujq- 
The imaginary part of the retarded local self-energy, 
Tiioc, then possesses a threshold ener gy, ri+ = Aq + uq, 
with ImS/oc l/v/(|f^| - wo)2 - Al for \n\ >nf, and 
ImEioc = for < $7+. As a result, ReEjoc exhibits 
a square-root singularity for |f2| —^ f2+ — 0^, which in 
the denominator of S leads to a suppression, i.e., dip, in 
the DOS. This dip is therefore similar in nature to that 
observed in angle-resolved photoemission experiments in 
the HTSC Finally, we find that the amplitude of the 
bound state peaks exhibits spatial oscillations described 
by cos{2kFr)e-''/yr, where I = kp/m^Al - Vll 

For non-zero temperatures the u) ~ +ujq branch of 
the bosonic mode becomes populated, opening a new 
scattering channel and yielding significant changes in 
the DOS, as shown in Fig. 133. We find that the spe- 
cific temperature evolution of the bound state peaks de- 
pends on whether f2f, at T = is smaller or larger than 
= Aq — LOf). This distinction arises since I]ij(fc,r2) 
now not only possesses poles at \Q\ = + ujq, with 
weight 1 + nB{uJo) (for T, wq ^ ^o), but also poles at 
|ri| = Ek — loq with weight nsii^o)- Hence, ImEjoc 7^ 
already for |f2| > fij, in contrast to the case at T = 0. 
Thus, a bound state peak with fi" < fif, < Ag at T = 
becomes damped with increasing temperature, as shown 
in Fig. I^b for T — O.ISAq. At the same time, a new 
peak appears in the DOS at \n\ = fij, since ImEjoc now 
also possesses square root singularities at ±^lj. With 
increasing temperature the damping of the bound state 
peak grows (T = 0.4Ao) and its amplitude consequently 
decreases, while the newly emerged peak moves towards 
lower energies. In contrast, if < i^^^ shown) 
the bound state peak simply moves towards lower ener- 
gies with increasing temperature, in agreement with our 
analytical results in Eq.Q. This downshift is expected 
since the population of the mode's branches grows with 
increasing temperature, thus leading to stronger scatter- 
ing. Thus, varying T or wq has an opposite effect on the 
energies of the bound state peaks in the DOS. 

We next consider a localized mode in a _y2 —-wave 
superconductor with SC gap Ak = Ao(cosfc2, — cosky)/2 
and Aq = 25 meV. To simplify the discussion we consider 
a particle-hole symmetric dispersion, ek = — 2t(cos k^ + 
cos ky) , with t = 300 meV. A band dispersion with strong 
particle-hole asymmetry, characteristic of the HTSC, 
yields qualitatively similar results • In contrast to the 
s-wave superconductor, a localized mode induces a reso- 
nance state even with a spin-independent coupling, ga = 
g. These resonances are a hallmark of the strong coupling 
theory, and are not observed in the weak scattering limit 
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FIG. 3: DOS at r = for r = 0+ and g = 1.2 eV. (a) T ^ 
and several values of uiq. Inset: arrows indicate the "dip" in 
the DOS. (b) loo = 10 meV and several T. (c) ujq = 20 meV 
and several T. 



[lot. In Fig. 121 we present the DOS at r = as a func- 
tion of cjo and temperature (F = 0+). For T, tJo = (see 
Fig. the energy of the resonance peaks, fir, is given 
by the solution of = ±27rAo/ {5 log [|f7r|/(4Ao)]}, 
similar to the energies of resonance states induced by 
a non-magnetic impurity Jsl (this result was also inde- 
pendently obtained by Si [l^l)- For > wq ^ 0, 
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(8) 



where K is the complete elliptic integral of the first kind, 
while ImE/oc = for \n\ < ujq. Thus there exist a 
real gap in ImSjoc, in contrast to ImG/oc- Performing 
a Kramers-Kronig transform, we find for f2/u;o 1, 
ReSfoc ~ —An/ojQ where A ~ dxIm'Sioc(x)/x'^. 
In the limit 17,. <C Ao, the resulting equation for Qr is 
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A log 



4Ao 
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(9) 



As expected, increasing luq leads to a shift of f2r to higher 
energies, in agreement with our numerical results shown 
in Fig. Moreover, for wo — > 00 we again recover the 
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FIG. 4: DOS at r = for T = and several values of T for a 
d^2_y2-wave and s-wave superconductor (see inset, the curve 
for r = 4Ao has been multiplied by 0.1). 

unperturbed DOS (see, e.g., ujq — 200 meV). The loga- 
rithmic divergence in Ejo^ at \Q\ — Aq+ujq leads to a dip 
in the DOS (see arrows in the inset of Fig. A much 
weaker dip emerges in the weak scattering limit ^3 • 

The specific temperature evolution of the resonance 
peaks again depends on the relative order of jri^ | at T 
and fl~ = A — loq, similar to the s-wave case. For \D,r \ < 
fl~ (see Fig. Et) the resonance peak is shifted to lower 
energies with increasing temperature. In addition, the 
opening of a second scattering channel for T ^ leads to 
a logarithmic divergence in Ejoc, and consequently a dip 
in the DOS at \n\ il^ ~ 15 meV, as indicated by the 
arrows. In contrast, for O J < \i}r\ < (see Fig. IS];), the 
same opening of a second scattering channel leads to an 
increase in ImE;oc for |r2| > wq, which in turn shifts the 
peaks to higher energies and suppresses their amplitude. 
The smaller peaks at \U\ = U~ = 5 meV (see arrows) 
arise from the opening of the second scattering channel, 
and the resulting logarithmic divergence of ImE/oc at n~ . 
Finally, we find that due to the momentum dependence 
of the SC gap, the spatial dependence of the resonance 
peaks (not shown) [ll| is qualitatively similar to that 
observed experimentally near impurities in the HTSC . 

To study the effects of a mode with finite lifetime, F 

0, we present in Fig.^ the DOS in a (i^2_y2-wave and s- 
wave (see inset) superconductor at T = and for several 
values of F. In both superconductors the peaks shift to 
higher energies with increasing F. Note that a non-zero F 
does not directly change the lifetime of the induced state, 

1. e., the width of the peaks. This is particularly evident 
for the s-wave case (see inset). A non-zero F lowers the 
threshold frequency for ImEjoc 7^ from ri+ = Aq + wq 
(for F = 0+) to fie = Aq. However, for < Aq, we still 
have ImEjoc = and the width of the bound or resonance 
peaks remains unaffected by F 7^ 0. 

Finally, we consider the lowest order vertex correc- 
tion, shown in Fig. ^p. For a quadratic normal state 



dispersion, the general form of this correction is Sg/g = 

— (^) F{^rn, Vn, i^o, Aq) whcrc i7„i, are the bosonic 
and fermionic Matsubara frequencies, respectively, and 
-F is a bounded function of 0(1). In particular, for an s- 
wave superconductor, we have F = -t-1 for Q.^, Vm^o 3> 
Aq, and F — —1 for flm^Vm^o ^ Aq, implying that 
the bare vertex is enhanced (suppressed) for frequen- 
cies smaller (larger) than the SC gap. In contrast, in 
a dx2-y2-wa,Ye superconductor we obtain ^ for 
rijrnh'nj^n *C Aq, and vertex corrections become irrele- 
vant. Thus, vertex corrections can in general be neglected 
for gm/A < 1, which applies to all cases considered above. 

In conclusion, we consider the effects of a localized 
bosonic mode on the electronic structure in s- and 
dx2_y2-weLYe superconductors. The mode acts as a new 
type of "defect" , leading to the emergence of bound or 
resonance states. We identify several characteristic fea- 
tures, such as a "dip" in the DOS and the temperature 
dependence of the mode induced peaks, that qualitatively 
distinguishes the effects of a localized mode from those of 
static impurities. This result provides further insight into 
the nature of impurities observed in STM experiments. 
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